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Abstract 



A manifestly invariant renormalization scheme of N = 1 nonabelian super- 
D ' symmetric gauge theories is proposed. 

43 



1 Introduction 



Construction of a manifestly supersymmetric and gauge invariant renormalization 
procedure is a nontrivial and in some cases not yet solved problem. Certainly, it is 
preferrable to have a calculation scheme which preserves supersymmetry at any inter- 
mediate stage. However, most popular gauge invariant regularizations like dimensional 
regularization [1], dimensional reduction [2] or lattice regularization break supersym- 
metry [3]. Higher covariant derivative regularization [4, 5] can be easily used in Abelian 
theories [6, 7, 8]. In principle, higher derivative regularization is also applicable to non- 
Abelian theories [9] , but the calculations of quantum corrections are quite involved. In 
this case it can be convenient to use non gauge invariant regularization. For example, 
one can regularize the theory by adding higher derivative term with ordinary deriva- 
tives, instead of the covariant ones for the Yang-Mills field which simplifies the calcu- 
lations considerably. Such regularization preserves supersymmetry, but breaks gauge 
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invariance, which should be restored by a proper choice of a renormalization scheme. 
For example, in the framework of algebraic renormalization it is done by tuning finite 
counterterms to provide relevant Slavnov- Taylor Identities (STI) for the renormalized 
Green functions. This method was applied successfully to SUSY gauge theories in the 
papers [10, 11], where the invariant renormalizability of N = 1 non-Abelian SUSY 
gauge models was proven in the framework of algebraic renormalization. However a 
practical implementation of the algebraic renormalization is rather cumbersome as the 
procedure requires a tuning of a large number of (noninvariant) counterterms. 

Recently a new method of invariant renormalization was proposed [12, 13], which 
provides automatically the renormalized Green functions satisfying STI. This method 
was generalized to sypersymmetric electrodynamics in [14] . However, in order to apply 
it to non-Abelian supersymmetric theories it is necessary to solve some problems, 
caused by a more complicated structure of STI. This is done in the present paper. 

The paper is organized as follows: 

In Section 2 we introduce the notations and remind some information about super- 
symmetric Yang-Mills theory. The universal invariant renomalization scheme for the 
model is constructed in Section 3. This scheme is illustrated in Section 4 by calcu- 
lation of the one-loop /3-function with the simplified (noninvariant) version of higher 
derivative regularization. The results are discussed in the Conclusion. 

2 N = 1 supersymmetric Yang-Mills theory. 

N = 1 supersymmetric Yang-Mills theory may be described by the following action: 

i 

S = ^Retr J d A x d 2 9 W a C ab W b + 

+ 16^ tr J d 4 xd 4 e(^{c + -c) ^5V + aB + B - iB + D 2 V -iBD 2 V^j, (1) 

which is invariant under the BRST transformations: 



5V = e 



-VctgV (c - c + ) + -V(c + c + ) 



2 

5B = 0; 5B + = 0; 

5c = ieB; 5c + = —ieB + ; 

5c = -ec 2 ; 5c + = -e(c + ) 2 . (2) 



1 In our notations the metric tensor in the Minkowski space-time has the diagonal elements (1, -1, 
-1,-1). 
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Here we use the following notations: 

V = -ieV a t a , ti(t a t b ) = ^5 ab , (3) 

where t a are hermitian generators, so that V a are real scalar superfields and for any 
function f(V) = c + c{V + c 2 V 2 + ... we define 

f(V)A = c A + Cl [V, A] + c 2 [V, [V, A}]+.... (4) 
W a is a chiral spinor superfield, given by the equation 

W a = LD(l-K)D(e- 2iV (l + K)D a e 2iV ), (5) 
where D is a supersymmetric covariant derivative: 

0=^-^0 8^ (6) 

It is convenient to introduce the generating functional with additional sources G 
and g: 

Z = J dfi exp jiS + itr J d A xd A 6 (yj + G-^SV) + itr J d 4 xd 2 6 x 

x{j c c + j c c + g^5c)+itr J d 4 xd 2 6 (j+ c+ +j+c+ + g +^-5c+)^. (7) 

Starting from this generating functional we construct the generating functional for 
connected Green functions 

W = -i\nZ (8) 

and the effective action 

r = w - J d 4 xd A evj - J d 4 x d 2 e (j c c + ] c c) - J d 4 x d 2 e (j+ c + + j+c+) , (9) 

where the sources should be expressed via equations 

TA SW 5W _ 5W . SW , 5W 

Note, that the effective action will depend on the sources G and g as on parameters. 
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3 Invariant renormalization 



Performing in the generating functional (7) a substitution (2), it is easy to derive 
STI, which can be written as 
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tr / <Fx<FB 



+ 



5V5G a\ 5c ' <?c+ 

+tr f fotfe^+tr I d 4 xd 2 9^-^- = 0. (11) 
J be bg J bc + bg + 



Due to invariance of the generating functional under the transformations c — > c + 
e(x), where e(x) is an arbitrary anticommuting function, the following identities take 
place: 

In order to write these STI in terms of the effective action it is convenient to define a 
functional T by subtracting the gauge fixing term from T: 

f = T ^-tr f d 4 x d 4 9 V (D 2 D 2 + D 2 D 2 ) V. (13) 

This functional satisfies the following equations: 



tr f d A xd 4 9 ^-^- + tr [ d 4 x d 2 9 6 -f 6 -f+tr [ d 4 x d 2 9 f-^- = 0; (14) 
J bV b& J be bg J bc + bq + 

i ^ sr sr „ i „ bf sf 



5c 5g J 5c + 5g^ 



^ d Jg + T, =0 - -3^ d Jg + JF = °- (15) 

Note, that up to now we did not set any field equal to 0. Moreover, identity (14) is 
not linear. So, to treat these identities the following algorithm should be used: 

1. In order to find relations between different Green functions it is necessary to 
differentiate identities (14) and (15) with respect to the arguments and then set all 
fields equal to 0. 

2. Identity (14) should be expanded over the Plank constant h, that corresponds 
to the loop expansion. Let us present the effective action as follows: 

f = S + hAT w + h 2 AT i2) + . . . (16) 

Then in the lowest order this identity is satisfied due to BRST-invariance of the classical 
action. At the arbitrary order one has 

n-l 

Ar (n) • s + s ■ Ar (n) = -J2 Ar (m_n) • Ar (m) (17) 

m=l 



4 



where • denotes an operator, which appears after differentiation of (14): 



A - B ^r xde [wsG + wsG) + 

+tr /^ A + __j + tr | Me + -— j = o. (is) 

The identity (17) holds if a regularization preserves gauge invariance. If a noninvariant 
regularization is used the invariance may be restored by adding local counterterms to 
the proper vertices. This is the procedure adopted in the algebraic renormalization 
scheme, where the main problem is to find the renormalized action S ren , providing 
the effective action which satisfies ST identities. As the explicit determination of all 
noninvariant counterterms is a difficult problem, we adopt another approach, which 
allows to write directly the gauge invariant renormalized effective action by solving ST 
identities. 

Our goal is to get the renormalized effective action T R , which satisfies STI by 
construction: 

f R = S ren + ftA 7 W + h 2 A^ + .... (19) 

Here A7W are some functions, which are constructed from Ar^ so that the renor- 
malized effective action satisfies renormalized STI. The renormalized action is given 
by 

S ren = S- hA^l - h 2 A^l + ..., (20) 

so that 

f R = S + hAT ( R V> + h 2 AT ( R Y> + ..., (21) 

where 

Alf = A7 W - A7& (22) 

If the theory is anomaly free, it is always possible to find such counterterms A^ v , 
that 

A7 W -A 7 « =Ar«-Ar«, (23) 

in the limit, when a regularization is removed. Here ArW and Ar^'„ correspond 
to A7W and A7^ in case of invariant regularization. Below we shall show that by 
solving STI one may get an expression for Ar^ without finding explicitely noninvariant 
counterterms A^ v . 
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We shall use the eq.(17) to determine the renormalized Green function of order n 
assuming that according to -R-operation the subtractions in all divergent subgraphs, 
corresponding to T^ m \ m < n are done in agreement with the eq.(17). 

We shall work with the equation (17) assuming that it is differentiated with respect 
to the fields, and then all the fields are put equal to zero. In this way we get a relation 
which expresses an arbitrary Green function with at least one chiral external gauge 
line in terms of other Green functions which are of the same order in h but have one 
external gauge line less and lower order Green functions. Symbolically 

D 2 x T n {x...) = R(T m ) (24) 

where R(T m ) denotes a product of lower order Green functions and the Green func- 
tions of the same order but having one external gauge line less. We assume that all 
Green functions of order m < n are renormalized in agreement with the eq.(17) and 
accordingly the divergent subgraphs in the function r n are subtracted following the 
R-operation. Then the eq.(24) fixes the overall subtraction in the function T n so that 
the renormalized function satisfies STL Indeed let us present the function T n in the 
form 

T n (9,p) = Y / B l (9,p)F l (p) (25) 

i 

where the functions Bi are polynomials in 9,p and Fj depend only on p. The renormal- 
ization is achieved by subtracting a polynomial of Fi(p). We will choose the polynomials 
Bi so that their chiral parts 

Q l = D 2 x B, t (26) 

are either linear independent, or equal to zero. 

The r.h.s. of eq.(24) may be decomposed over the polynomials Qi 

R(r m ) = -£Qi(°,p) R i(p) (27) 

i 

In practice to find the functions Ri explicitely one has to perform a multiplication of 
the coefficients functions of the Green functions T m , entering the r.h.s. of eq.(24) using 
the supersymmetry algebra. 

According to our assumption all functions at the r.h.s. of eq.(24) are finite. Using 
linear independence of coefficient functions Qi, we find that the renormalized function 
T„ may be written as follows 

It = J2B i (9 : p)R i (p) + J2B j (9,p)F j (p) (28) 

i 3 

where the second sum runs over j for which Qj = 0. It completes the renormalization. 
Such renormalization can be made in the following sequence: 

At the first step we should calculate one-loop correction to the two-point Green 
functions of the ghost fields: 
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-Ar 



a). 



aFST"" w Ar ™ (29) 

STI do not impose any restrictions on the renormalization of these functions. There- 
fore, they can be renormalized in arbitrary way say by subtraction at some fixed 
normalization point /x c . At the next step we calculate functions 



52 Art 1 ). 



52 ArW 



(30) 



For this purpose it is necessary to use STI, which is obtained by differentiating equation 
(15) over c. Then we proceed to the renormalization of the function 

a 3 Arw 

s^sEsv' (31) 

using the corresponding STI. Having constructed functions (30) and (31) it is possible 
to renormalize the function 

£ 3 Arw 

5c+ 5c5V ^ 
In each case the STI, which are needed for the renormalization, are obtained by dif- 
ferentiating of (14) or (15). Then we differentiate all STI over V and proceed to the 
renormalization of functions with larger powers of V, starting from 

5c5G5V (33) 
Renormalization of the two-point Green function for gauge field is made similar to 
the QED case [14]: If we introduce the function II 



(2tt) 4 5 4 (( Pi )m + ( P2 ) M ) u[(e xl , Pl ), (e X2 , - Pl ) 



J d 4 xid A X2 



5 2 T 



SV X1 5V X2 



<t>,<t>,v=o 



exp (i( Pl )^ + i(p 2 )^2 ) , ( 34 ) 



then 



n 



(e x , P ), (e y , - P ) = f 1 ( p ) p 2 n 1/2 5\e x - 9 y ) + F 2 (p) 5\e x - e y ), 



(35) 



where 



n V2 = -—D a D 2 C ab D b . 



lQd 2 ' 



In this case 



(36) 
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B 1 (9,p)=p 2 U 1/2 5\9 x - 6 y ); B 2 (9,p) = 8\9 X - y ) 



(37) 



The renormalized Green functions are obtained by subtracting from F{{p) some 
polynomials chosen to provide STI for the function IF 

TT[6,p] = J2B l (6,p)(F l (p) - P(p)), (38) 

i 

For the two-point Green function substitution of (35) into the STI 



D 



y 5V z SV y 



= 0, 



(39) 



which is obtained by differentiating the linearized STI (14) over c + and V, gives the 
following equation: 



(F 2 (p) - P 2 (p)) D 2 X 5\6 X - 6 y ) = 0. 



(40) 



Therefore, 



P 2 (p) = F 2 (p), (41) 

while the function Pi can not be defined from STI and corresponds to a gauge invariant 
counterterm. It is convenient to choose 

Pi(p) = Fi(/U (42) 

where normalization point. Then the renormalized two-point Green function 

can be written as 



U r [(6 x ,p), (6 y , -p)] = (F 1 (p) - F 1 (^))p 2 Il 1/2 5 A (9 x - 6 y ). 
This Green function satisfies the equation 



(43) 



DlW 



e ,p),(e y ,-p)\ =0, (44) 

which is a supersymmetric generalization of transversality condition in the usual quan- 
tum electrodynamics. 

Renormalization of the three-point function for the gauge field can be made accord- 
ing to the above algorithm using the corresponding STI. 

The procedure described above may be interpreted as finding local counterterms, 
which can be written as 

1 r / oo oo \ 

AS=— tr J d 4 xd 4 9(c + -c)UY / a n V n (c-c + )+ ]T b n V n (c + c + ) + 

\ n=l n=0,7^1 / 

+tr [ d 4 xd 4 6G(if^f n V n (c-c + ) + f^h n V n {c + c + ) \ + tr f d A xd A 0F[V}. (45) 

•* V n=0 n=0 ) J 
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After adding these counterterms the renormalized effective action Ar^ will satisfy 
renormalized STL However, the above algorithm allows to avoid writing and tuning of 
all terms in (45), that may cause considerable technical difficulties. 



4 One-loop renormalization with the simplest PV- 
regularization 

In order to illustrate the scheme, presented above, we consider N — 1 supersym- 
metric Yang-Mills theory and calculate its /^-function using regularization by (usual) 
higher derivatives, complemented by Pauli-Villars fields for regularization of one-loop 
divergences. 

The regularization is introduced by adding the higher derivative term, so that the 
regularized action is equal to 

1 r ( 8 2n \ 

S reg = ^Retr J d A x d 2 6 W a C ab il + — \W b + 

+Y^tr J d 4 xd 4 e(^c + -c)-^5V + aB+B - iB + D 2 V -iBD 2 V^j, (46) 

However, the higher derivative term does not regularize one-loop divergences. There- 
fore, it is necessary to add in the generating functional PV-fields, that can be made by 
the following way: 

Z = J rf^ [[ ( det M,))"' [I ( det mi )) ft x 

i i 

x exp ^iS reg + itrj d 4 x d 4 6 (VJ + G ^5V) +itrj d A x d 2 6 x 

x(jcc + ] c c + g-^5c) +ttr J d 4 xd 2 6 (j+ c+ +j+c+ + g+ j^c+) j. (47) 

where 



(det PV(V,M)) 1 = 

= / DWe^Uj d'xd%d'yd%W:-^- b W;- l -M 2 J d*xd*9W 2 \, 

\ x y J 

(detpi>(V, to)) = 

= J DC DC exp | ^tr J d 4 x d A o((C + - C) ^5V(V,C,C + ) + 
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+ !J tr J d 4 xd 2 9CC + ^tr J d A xd 2 6C + C + ^ (48) 
and the coefficients cej and fa satisfy the following equations: 

5>i = i ; E^ M ' = °; 

Eft = 1 ; EK = °- ( 49 ) 

Here W a , C and (7 are PV fields, Mj are masses of PV-fields, corresponding to V 
superfield and rrii are masses of PV-fields, corresponding to ghosts. We will assume, 
that all Mj and are proportional to the parameter A in the higher derivative term. 

In order to find two-point Green function for the gauge superfield it is necessary to 
calculate Feynman diagrams, presented at Figure 1. The result appeared to be 



= -*e 2 Cl / ^- V a (-p) d*n 1/2 v«( P ) x 



x ■ 



I(i + (-i)V7A 2 ") / 



d A k 



(2vr) 4 \k 2 (k + p) 2 (l + (-l) n k 2n /A 2n ) 



y^^fc 7 w r ] + 

k (k 2 - M 2 ) Uk + p) 2 - M 2 ) (1 + (-l) n k 2n /A 2n ) J 



1 t d 4 k ( 1 v r d 4 k 1 \1 

+ ll j2nY\k 2 {k + p) 2 \ ak l (2nY {k 2 _ M 2 ) JJk + p) 2 -M 2 ))} + 

d 4 p 



+e% J j^- 4 V a (-p)V a (p)x 



(2k) 

j^iaj f d 4 k 1 ^ iocj r d A k 2(k + p) 2 + p 2 

x l 2^ ~o~ / io^\i (U2 _ n/r2\ ~ 2^ ~q~ / ToT/u 777 7777777 7^ 7777 ^ 



+ E 



j 

i a 



3 

" 4 

j 



3 J (2nY{k 2 -M 2 ) ^ 8 J (2n)^ k 2_ M 2 )((k + p) 2 -M 2 ) 
d 4 k M 2 {l + (-ir(k + p) 2n /A 2n ) 



4 J (2VT) 4 (Jfe2 _ M 2) (l + (_l)np«/A2n) ((jfe + p )2 _ M ^ 



10 / f<W\4 fl.2 _ ™2\ o / /r. \4/,o , _\o o\ I ^ ' 



12 J (2tt) 4 (P - mj) Y 8 7 (2tt) 4 _ m 2) ^ + p)2 _ m ^ 

Substituting this expression to equation (34) we obtain the function IT, which is equal 
to (35) with 



^) = ^{( 1 + (- 1 )V7A 2n )/(0i( 



k 2 (k + p) 2 (i + (-i) n k 2n /k 2n ) 
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^ (k 2 - Mf) ((Jfe + p) 2 - Mf) (l + (-l) n k 2n /A 2 ^ ' 



+ 



i /, ■ — — ■ i ! - ) i i /.. • — : - ■ ,' i ; ■■ - -- ■ t ' — i i . — i— i — . )"/■." — ■ 1 / \ — ' ' 

rf 4 /fc / 1 ^ /■ d A k 1 



I r d*k 1 v /• 



(2tt) 4 VA: 2 (A:+p) 2 V * •/ (2tt) 4 _ M 2 ) ((A; +p) 2 - M 2 
d 4 fc 1 



ia 



ia 



, ■ . . , d'k 2(k + p) 2 + /r 

' ~ ' (2tt) 4 (P_m 2 )((A; + j9) 2 -M 2 ) + 



^ M 2 (l + (-l)"(A; + p) 2 7A 2 «) 



+ ^ 4 i (2vr) 4 (A;2_M 2 )(l + (-l)™A; 2 ™/A 2 ™)((A; + p) 2 -M 2 ) 

i/^ /■ tf*fc 1 r d A k p 2 1 

^ 12 i (2tt) 4 (A; 2 - m 2 ) + ^ 8 7 (2tt) 4 (p _ m?) ( (jfc + p y _ m ^ } ^ > 

In order to calculate the integrals in this expression, first it is necessary to perform 
the Wick rotation. Then it is easy to see, that 

3 A 

Flip) = e2 Ci( m — h finite terms) (52) 

16 v p ' 

It is easy to verify, that the function F 2 is finite and proportional to p 2 . After taking 
a limit A — >■ oo (and therefore, Mj — >■ oo and mj — >■ oo) we obtain, that 



n 



3 A 
%,P), {0 y , -P)\ = —r^e 2 c 1 p 2 U 1/2 5^(e x - 0j,)(ln- + finite terms) + 



167T 2 ^ ±/Z V ^ ~»>\-p 

+e 2 Cl 5\6 x - 6 y )p 2 C( aj , (3j, Mj/A, mj /A) j. (53) 
where the constant C is given by 

p 2 C = lim F 2 (p) (54) 

A^oo 

and in the general case is not 0, that can be easily verified. Then the renormalized 
function IT, obtained according to prescription (43), is 



n 



{0 x ,p),{9 y ,-p) = -T7 —e 2 c 1 p 2 U l / 2 5' l (d x - 0„)( In- + finite terms) (55) 



J- 2 e 2 c lP 2 U 1/2 5\9 x -9 y )(\n^ 

and corresponds to the following renormalized contribution to the effective action: 
3e 2 f d A p 



= 32^ Cl / (2^ ^"^(p) ( ln ~ v + finite terms ) • ( 56 ) 
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One-loop correction to the ghost propagator corresponding to the first diagram, 
presented in Fig. 2 is found to be zero. We also verified, that the one-loop ghost-gluon 
vertex, corresponding to the other diagrams in this figure, is finite. Therefore, the 
one-loop /5-function 

^^i&M)' (57) 

calculated in the noninvariant regularization, described above is 

(3(a) = + 0(a 3 ) (58) 

and agrees with the calculations made by DRED. However, it is interesting to find the 
two-loop contribution to the /^-function (57) in the considered regularization, because 
in SUSY QED [15, 16] such contribution is different from the DRED result and is equal 
to zero. 



5 Conclusion. 



In this paper we presented a renormalization procedure for SUSY Yang-Mills theo- 
ries, which guarantees BRST-invariance of the renormalized theory for any intermedi- 
ate regularization. SUSY Slavnov- Taylor identities are incorporated into subtractions, 
which allows to avoid the appearance of noninvariant counterterms. The proposed 
procedure is illustrated by the calculation of the one-loop /3-function for the N — 1 
SUSY Yang-Mills theory with the noninvariant version of higher derivative regulariza- 
tion. However it would be interesting to calculate two-loop quantum corrections and 
compare them with the correspondent result in SUSY QED [15, 16]. 
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